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Abstract 



We analyse the anisotropy of homogeneous turbulence in an electrically conducting fluid 
submitted to a uniform magnetic field, for low magnetic Reynolds number, in the quasi- 
static approximation. We interpret disagreeing previous predictions between linearized theory 
~~^ and simulations: in the linear limit, the kinetic energy of transverse velocity components, 

__^ normal to the magnetic field, decays faster than the kinetic energy of the axial component, 

r^ along the magnetic field (Moffatt (1967)); whereas many numerical studies predict a final 

^^ state characterised by dominant energy of transverse velocity components. We investigate the 

■ O corresponding nonlinear phenomenon using Direct Numerical Simulations of freely-decaying 

turbulence, and a two-point statistical spectral closure based on the Eddy Damped Quasi- 
Normal Markovian model. The transition from the three-dimensional turbulent flow to a 
"two-and-a-half-dimensional" flow (Montgomery & Turner (1982)) is a result of the combined 
^ effects of short-time linear Joule dissipation and longer time nonlinear creation of polarisation 

anisotropy. It is this combination of linear and nonlinear effects which explains the disagree- 



^ ment between predictions from linearized theory and results from numerical simulations. The 



transition is characterized by the elongation of turbulent structures along the applied magnetic 
field, and by the strong anisotropy of directional two-point correlation spectra, in agreement 
with experimental evidence. Inertial equatorial transfers in both DNS and the model are 
presented to describe in detail the most important equilibrium dynamics. Spectral scalings 

^~H are maintained in high Reynolds number turbulence attainable only with the EDQNM model, 

^ which also provides simplified modelling of the asymptotic state of quasi-static MHD turbu- 

^fi lence. 

m 

■^ Magnetohydrodynamics, Quasi-static hypothesis. Homogeneous turbulence. Direct Numerical 

Cn Simulations, EDQNM 

^H 1 Introduction 

^ In most geophysical and astrophysical flows, turbulence is affected by forces that distort signifi- 

• 'H cantly some of its scales in an anisotropic manner, such as the Coriolis force in rotating flows or the 

rS Lorcntz force arising from the presence of an external magnetic field in a conducting fluid. This 

3 specific turbulent dynamics forced by an imposed magnetic field is found in liquid metal flows, be 

they of industrial, geophysical nature — the melted iron core of the earth — or of academic interest 
in the laboratory, such as the experiment by Alemany et al. (1979) in liquid mercury. Recent lab- 
oratory experiments on the dynamics of conducting fluids use sodium or gallium; liquid sodium is 
also used in industrial configurations, for instance in the French fast breeder reactor Superphenix. 
Generally, the motion of turbulent liquid metals is governed by magnetohydrodynamics (MHD): 
the induction equation for the fiuctuating magnetic field is added to the Navier-Stokes equations, 
which are in turn modified by the Lorentz force, representing the feedback from the magnetic 
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field. In the presence of an external magnetic field, such MHD coupling results in new dissipative 
terms, of ohmic nature, and selectively damped waves, the Alfven waves (Moffatt (1967)). In 
cases involving liquid metal, the magnetic difFusivity in the induction equation is larger than the 
molecular viscosity in the Navier-Stokes equations, i.e. the magnetic Prandtl number is small 
compared to one. The magnetic diffusivity is so large with respect to the kinematic difFusivity — 
with a magnetic Prandtl number less than 10^^ in the Earth's iron core, of order 0.9 x 10^"'' in liquid 
sodium or 1.4 x 10~^ in mercury — that it is consistent to consider the flow at very high Reynolds 
number and at low magnetic Reynolds number. In the following simulations, the magnetic Prandtl 
number is set to Ptm — 3.1 x 10^'*. 

As discussed in section 2, if the magnetic Reynolds number is small enough, the linear regime no 
longer admits Alfven waves solutions, and the effect of the Lorentz force reduces to an anisotropic 
ohmic (or Joule) dissipation term. In this regime, called the quasi-static approximation (QS MHD), 
the induction equation is simple enough to be solved explicitly and to yield a closed expression of 
the Lorentz force in terms of the velocity. The specificity of the quasi-static limit can be discussed 
both in terms of timescales and anisotropy. Unlike more general MHD turbulent flows, in which 
nonlinear and Alfven timescales may be in competition and yield length scale dependent levels of 
anisotropy (see e.g. Zhou & Matthaeus (2005); Zhou (2010)), in QS MHD the magnetic diffusivity 
is too large to enable Alfven waves. The only relevant timescales concern the modified Navier- 
Stokes equations, with a linear nondimensional timescale tj/Bq resulting from ohmic dissipation 
(Bq is the external magnetic field, scaled as velocity, and rj the magnetic diffusivity), and the 
nonlinear timescale Io/uq (uq is the rms velocity and Iq the length scale related to a turnover 
time). Strong anisotropy is first induced by the ohmic dissipation term over the linear timescale. 

Quasi-static MHD turbulence was investigated experimentally by Alemany et al. (1979) and 
Caperan & Alemany (1985). In these studies, turbulence was generated by towing a grid through 
a cylindrical tank full of mercury, with an external magnetic field generated by a coil. Measure- 
ments include Reynolds stress components, an integral length scale in the axial direction and 
one-dimensional spectrum of transverse energy with respect to the axial wavenumber. A clear 
transition from a three-dimensional state, with conventional Kolmogorov spectrum, to a quasi- 
two-dimensional state, with fciT spectrum, was evidenced. The first phase of this 3D-2D transition 
was studied using axisymmetric Lin equations with an Eddy Damping Quasi-Normal Markovian 
(EDQNM) closure model by Cambon (1990), and the scenario of a two-dimensionalization in two 
steps was proposed. This scenario was recently confirmed by DNS in Favier et al. (2010) and one of 
the goals of the present paper is to go beyond the numerical approach of Cambon (1990) using both 
anisotropic EDQNM and direct numerical simulations (DNS). The "eddy-damping" rate appearing 
in the EDQNM closure for general MHD turbulence should in principle be modified to account for 
the combination of sweeping and straining mechanisms, thus allowing for the possibility of either 
Kolmogorov inertial scaling {k~^'^ kinetic energy spectra, isotropizcd, i.e. spherically integrated) 
or Iroshnikov-Kraichnan scaling (fc~^/^) (see e.g. Zhou et al. (2004)). The QS MHD approxima- 
tion, without Alfvenic propagation, allows to anchor the model within the classical hydrodynamic 
turbulence context, thus keeping the original damping consistent with Kolmogorov scaling. 

Other numerical approaches in the same context are given by Schumann (1976), and Knaepen 
et al. (2004) with application to anisotropic modelling. A survey is offered by Knaepen & Moreau 
(2008), in which the change of anisotropic structure for the Reynolds stress tensor, from purely 
linear to nonlinear dynamics, is presented as an open problem. We think that this problem can be 
elucidated by the scenario of 3D-2D transition in two steps (Cambon (1990); Favier et al. (2010)) 
which is fully described hereafter. 

Both spectral theory and DNS were applied by Ishida & Kaneda (2007) to the dynamical and 
structural study of the small scales anisotropy of QS MHD turbulence, while a recent approach by 
Okamoto et al. (2010) focused on the infrared limit, i.e. at very large scales. In the latter work, 
assuming the existence of a Loitsyanski-like invariant, decay laws for typical integral lengthscales 
and Reynolds stress components are proposed and compared to DNS results. The dynamics of 
integral length scales was shown to be crucial in rotating turbulence which bears strong analogies 
with MHD turbulence. For instance, the linear growth rate of the integral length scale related 
to transverse velocity components and axial separation, denoted £n in Okamoto et al. (2010), 

and L\^ — L22 here, was clearly related to the role of nonlinear transfer terms (Cambon & 
Jacquin (1989); Jacquin et al. (1990); Cambon et al. (1997)). This result was recently recovered 



by Staplehurst et al. (2008) with a different interpretation, although we beheve that the use of 
axisymmetric Lin equations — equations for two-point velocity correlation spectra (von Karman 
& Lin (1949)) — , in which linear and nonlinear terms are exactly separated, is essential to the 
understanding. Accordingly, our theoretical approach is based on an anisotropic spectral formalism 
with generalized Lin equations instead of on a formalism based on the Karman-Howarth equation, 
rather used by Okamoto et al. (2010), but bridges between the two approaches will be discussed 
in the following. 

One of the most challenging aspects of quasi-static MHD turbulence, from a numerical point 
of view, is the rapid increase of the velocity correlation lengths in the direction of the imposed 
magnetic field. In that case, the results from classical pseudo-spectral methods with periodic 
boundary condition are often questionable, as the characteristic scale of the turbulent motion is 
no longer small compared to the numerical box size. In this paper, we compare Direct Numerical 
Simulations (DNS) with a model based on EDQNM closures and confirm that neither the low 
Reynolds numbers considered in DNS nor the confinement due to periodic boundary conditions 
alter our understanding of the dynamics. Secondly, our goal is to propose a detailed study of 
the anisotropy of quasi-static MHD turbulence at low, moderate and high Reynolds numbers. 
As in Favier et al. (2010), the analogy with the asymptotic quasi-two-dimensional state, called 
"two-and-a-half-dimensional" flow, will also be discussed. 

The paper is organised as follows. The main parameters and governing equations are recalled in 
the following section. Spectral properties and EDQNM closures are discussed in section 3, and the 
numerical methods used in the paper are presented in section 4. Section 5 is devoted to the issue 
of confinement, both in DNS and EDQNM. Most of the results are gathered in section 6, where the 
statistical properties of quasi-static MHD turbulence are described, with an emphasis on anisotropy 
characterisation (section 6.2). Finally, the large Reynolds number behaviour is investigated in 
section 7.1, along with the analogy with quasi-two-dimensional turbulence in section 7.2. Details 
about EDQNM closed equations and linear predictions for the velocity correlation lengths are 
gathered in Appendices A and B. 

2 Governing equations and parameters 

We consider initially isotropic homogeneous turbulence in an incompressible conducting fluid, 
in which u^ — Uy — Uz, where u^^ Uy and u^ are the rms values of the velocity components. 
When the external magnetic field is applied, along z in the following, u^ will be called the axial 
component and Ux^ Uy the transverse components. The fluid is characterised by the kinematic 
viscosity v, density p and magnetic diffusivity 77 = (cr/io)^^; u is the electrical conductivity, /ig 
the magnetic permeability. These physical properties are assumed to be constant. The integral 
length scale is Iq, defined from the two-point velocity correlation tensor Rii{r) — {ui{xi)ui{xi + r)), 
as Iq = L Rii{r) / Rii(0)dr , (or equivalently from the kinetic energy spectrum). The Reynolds 
number and its magnetic counterpart are Re = (uqIq)/v ^ 1 and Rm ~ (""o^o)/?? ^ 1- The ratio 
between these two numbers defines the magnetic Prandtl number PrM = ^^/Vj which is very small 
in our study. The fiow is submitted to a uniform vertical magnetic field B scaled as Alfvcn speed 
as Bq — B/y/pfi(j. The ratio between the eddy turnover time Iq/uo £^nd the ohmic time tj/Bq 
is the magnetic interaction number TV = (B'^Iq) / {tjuq) . Within the quasi-static approximation, 
which implies that Rm tends to zero, but which is nonetheless approximately valid for all Rm < 1 
(Knaepen et al. (2004)), the Navier-Stokes equations become 

du „ 1„ „9 „o , 1 9^M 

—-+u-Vu = — Vp + vV^u + M^A-^—^ 1) 

at p oz-^ 

F 

where F is the rotational part of the Lorcntz force, A^^ is the inverse of the Laplacian operator, 
Mq = B^/rj and z the axial coordinate, along the direction of Sq. Compressible effects are not 
taken into account here, so that V • m = 0. 



3 Exact and model equations for two-point second-order 
statistics 

We obtain hereafter the equations for the spectral statistics of the second-order moment of the 
fluctuating velocity field u. The derivation is facilitated in two ways: first, by beginning with 
the Fourier coefficients of u before computing the second-order moments; second, by using a 
Helmholtz-like decomposition in order to derive all the algebra only in terms of the incompressible 
components, namelly the toroidal/poloidal decomposition. 

Equation (1) for the velocity is 3D-Fourier transformed, with Fourier coefficients denoted with 
'^~, and the pressure term is eliminated using incompressibility, introducing Kraichnan's projector 



^inin\^) — c\ 






^2 ; ' ■-« y^rn ^2 



(2) 



so that 



— +Vk'^ +MlcOS^e\ Ui{k,t) = Pirnnik)u^!^, (3) 

where k is the wave vector and 6 its orientation with respect to the z-axis. The unique new term 
reflecting the quasi-static MHD effect is algebraic, (Mq cos^ 8) Ui{k, t). 

For second-order velocity correlations, the most general information is given by the second-order 
spectral tensor Rij{k,t) which in the homogeneous case is given by 

{u*{p,t)u,{k,t)) ^ R,j{k,t)5\k - p) . (4) 

The 3D Dirac function expresses that only the Fourier velocity components at the same wave 
vector have non zero double correlation. Another expression is obtained by considering a discretized 
velocity field, as in DNS (thus turning the mathematical formalism of distributions and generalized 
integrals, applied in continuous space, to classical integrals applied to discretized functions). For 
the particular case of a cubic periodic domain of size L, this replaces the Dirac term in the above 
equation by a factor (L/(27r))^. 

The brackets in equation (4) denote statistical ensemble averaging: in DNS started with a single 
realization of the velocity field, statistical averaging is obtained by spatial averaging, assuming 
ergodicity and using the particular symmetries preserved here, namely axisymmetry. 

In the quasi-static MHD case under consideration, statistical symmetry is thus restricted to 
axisymmetry with mirror symmetry (the mean helicity is zero if initially zero), and the spectral 
tensor can be expressed in terms of toroidal and poloidal components of the velocity field in 
Fourier space. The two components are obtained using a polar-spherical frame of reference with 
base vectors e(^)(fc) and e'^^^(fc) {a.k.a. Craya-Herring frame of reference, see figure 1; Herring 
(1974)), as 

M(fe,f) ='u(i)(fe,t)e(i)(fc)+w(2)(fc^t)e(2)(fc) . (5) 

This decomposition automatically treats the velocity field as solenoidal, i.e. divergence free in 
physical space, through the algebraic orthogonality condition k ■ u{k) = 0. In addition, it allows to 
construct any related statistical correlation, with a minimal number of components, for arbitrary 
anisotropy. The decomposition (5) is general, although, since it relies on the arbitrary choice 
of a polar axis, it is especially well suited to axisymmetric configurations, in which the tensors' 
dependence reduces to the wavenumber k and its angle 9 to the axis. 

The expression for the two-point second-order spectral tensor is therefore 

i?,^.=$ieWeW+$2ef)ef , (6) 

in which all the tensors and vectors depend on k and 9 as the toroidal and poloidal energy tensors 

$^(fc,t) ==:$^(fc,cos0,i), ^'^{k,t) = ^'^{k, con 9, t) . (7) 

Considering the symmetries of the fiow, the most general decomposition in terms of energy density 
e, polarization Z and helicity % reduces to (Cambon & Jacquin (1989)) 

e(fc,cos6',i) = ^ ($^ + $2)^ Z(fc,cos0,t) = ;^($2-$i), n{k,t)^0. (8) 




Figure 1: Craya-Herring frame {e^^\ e^'^\ e'-^)) in Fourier space. In the general case, Fourier modes 
in the blue region contribute to E{k,9) (eq.(25)). However, if k is vertical, the sum of the two 
components u^^\k) and u(^''(fc) generates a vertically sheared horizontal flow (VSHF), and if k 
is horizontal, they correspond to transverse and axial components. Therefore, the polar modes 
{9 ~ 0) contribute to horizontal kinetic energy, whereas equatorial modes {9 ~ 7r/2) contribute to 
both axial (along e*^^)) and transverse (along e*^^^) kinetic energies. 

The polarization term Z is in general complex-valued and its imaginary part corresponds to a non 
zero cross-correlation between poloidal and toroidal velocity components. Here, Z is real-valued 
and both sets of statistical quantities, ^^ and $^, or e and Z are equivalent. 
It is straightforward to derive the following exact equations for e and Z: 

— + 2vk^ + 2Ml cos^ 6* j e(fc, cos 9, t) = T^^^ (fc, 6*, t) (9) 

— +2vk^ + 2Mlcos^ 9\z{k, cos 9, t) =T^^\k,9,t). (10) 

These equations are exact in the limit of homogeneous quasi-static MHD turbulence. They gener- 
alise the Lin equation, with the definition of cubic T^*^'^^ terms given in Cambon & Jacquin (1989), 
and recalled in appendix A. All the terms in these equations can be obtained in pseudo-spectral 
DNS, as in Favier et al. (2010), using summation of Fourier modes on rings, in contrast with the 
summation of Fourier modes on spherical shells as usual in the analysis of isotropic turbulence (see 
figure 1). However, the anisotropic (fc, 9) distribution of T^^^^^^ is more affected by lack of sampling 
and noise in DNS, especially at small k where /S.k/k is large. It is therefore worthwhile to develop 
a model based on equations (9) and (10) to evaluate the behaviour of the second- and third-order 
moments — energy and energy transfer spectra. The model may then provide smooth values for 
these quantities, to be quantitatively compared to DNS results. 

We will be using hereafter such a model, drawn from the anisotropic EDQNM closure the- 
ory, which has already been successfully applied to rotating or stably stratified turbulent flows, 
including a comparison with DNS (Cambon et al. (1997); Godeferd & Staquet (2003)). In the 
derivation of the model, the toroidal/poloidal decomposition proves useful and valuable for simpli- 
fying the expressions for triple velocity correlations, without using projection operators inherited 
from equation (3). Another simplification comes from the use of a slightly modified decomposition 
of velocity, u{k,t) = ^+(fc,i)7V(fc) + ^-{k,t)N*{k) analogous to (5), which brings out the helical 
modes ^± by projection onto N{k) ~ e^^)(fc) — ie(-'^)(fc) and N*{k) = N{—k). Helical modes are 
advantageous because they diagonalise the curl operator and allow a more compact decomposition 
of triple velocity correlations at three points (triadic terms), even in isotropic turbulence (see for 
example Waleffe (1992)). [In rotating turbulence, the helical modes are also the inertial waves 
modes, Cambon & Jacquin (1989); Waleffe (1993); Cambon et al. (1997); Belief et al. (2006)]. The 
starting point of the closure is the third-order spectral tensor S related to helical modes, defined 

by 

(6" (g, 06' (P, m {k, t)) - Sss's" {k, p, t)S^{k + p + q). (11) 



The generalised Eddy Damping Quasi-Normal (EDQN) technique is then applied to the equation 
that governs the third-order spectral tensor, 



- + u{e+p^ + q^) 



+ M^ {cos^ Ok + cos^ 9p + cos^ dq) Sss's"{k,p,t) = nss's"{k,p,t), (12) 



in which diss's" {k^p, t) represents the contribution of fourth-order velocity correlations. In order to 
obtain a closed set of equations, S2 is expressed in terms of sums of products of double correlations. 
This would be an exact evaluation of the fourth-order moments, were it applied to a Gaussian 
random variable (the 'QN' part). We apply a corrective term (the 'ED' part) due to the non- 
vanishing fourth-order cumulant, to account for the departure from Gaussianity of both third- 
order and fourth-order cumulants. We shall use the version of the model that has provided the 
best results in rotating or stably stratified turbulence. This EDQNM2 model, say, accounts for the 
anisotropic Joule dissipation in both the second-order moments equation and in the third-order 
moments one (12). When informative, the results of EDQNM2 will also be contrasted with those 
of the simpler EDQNMl model, which retains the Joule dissipation term only in the second-order 
moment equation, discarding it in equation (12). Contrasting both models allows to tell whether 
the main anisotropic mechanism is mostly linear or nonlinear. Additional information on the 
models is given in appendix A. 

4 Numerical methods 

To assess the validity of EDQNM closure in the context of quasi-static MHD turbulence and to 
obtain results at low and moderate Reynolds numbers, we perform Direct Numerical Simulations 
of equation (1) using a pseudo-spectral method implemented on a parallel computer. The velocity 
field is computed in a cubic box of side L with periodic boundary conditions using 512^ Fourier 
modes. [The conventional shorthand relationship L — 2tt for non-dimensional DNS is used here, 
except for the previous discussion after equation (4).] A spherical 2/3-truncation of Fourier modes 
is used to avoid aliasing and the time scheme is third-order Adams-Bashforth. The dissipative 
viscous plus ohmic terms are treated implicitly. 

The DNS results presented here are performed at higher resolution than those of Favier et al. 
(2010). An initially isotropic turbulent velocity field is created by a hydrodynamic simulation with 
large-scale forcing in order to reach a quasi-steady state. At the end of this prc-computation stage, 
the rms velocity is uq — 0.81 and the integral scale ?o = 0.25 yielding Re — uqIo/v ~ 333. The 
Reynolds number based on the Taylor microscale is R\ ~ 95. This rather low value, considering 
the resolution, is a consequence of our specific choice of a small initial integral scale Ig, in order 
to lift partially the numerical confinement constraint, discussed in section 5. The corresponding 
turbulent flow fleld is used as initial state for two different MHD simulations. In all of them 
Rm — 0.1 (hence A^ ~ 2), so that the quasi-static approximation is justified (Knaepen et al. 
(2004)). Two different amplitudes of the imposed magnetic field are chosen, which correspond to 
two values of the interaction parameter: TV = 1 and 5. For reference, we also compute the isotropic 
case, setting Bq — 0, from the same initial condition. The quasi-static MHD simulations are freely 
decaying to avoid spurious effects of a forcing scheme on the development of anisotropy. 

The two versions EDQNMl and EDQNM2 mentioned above are used. The EDQNM sim- 
ulations are initialised with the exact initial kinetic energy spectrum obtained from the DNS 
pre-computation. The EDQNM spectral space is discretized as follows: we use 64 values for the 
wave number fc, 32 for the polar angle 9, and 32 internal orientations for the angle defining the 
orientation of the plane of the triads. In contrast to DNS, the wave number discretization used 
here is logarithmic, thereby improving the sampling of the large scales with respect to DNS. The 
minimum and maximum wave numbers solved are kmin — 1 and fcmax = 512/3 as in DNS. If Bq = 0, 
EDQNMl and EDQNM2 are identical, and we also compute this particular case for comparison 
with isotropic DNS. Hereafter, DNS results are plotted with lines only, EDQNM results are plot- 
ted with lines and symbols (o for EDQNMl results, • for EDQNM2 results and A for isotropic 
EDQNM). 




t* = tuo/k 



Figure 2: Velocity correlation lengths versus dimensionless time t* = tuo/lo. The thick lines 
correspond to DNS without nonlinear interactions. The thin lines correspond to analytical linear 
predictions from Appendix B. The horizontal line presents the numerical limit of 27r. 

5 Confinement due to periodic boundary conditions 

This paragraph is specifically devoted to the problem of confinement in quasi-static MHD turbu- 
lence. As the anisotropic ohmic dissipation affects the fiow, the velocity field rapidly homogenizes 
in the direction of the imposed magnetic field. The velocity correlation lengths thus increase in 
the axial direction. However, due to periodic boundary conditions used in DNS, these correlation 
lengths are limited by the size L = 27r of the computational domain. To remove possible non phys- 
ical effects due to this confinement, we compute the initial velocity field with an integral length 
scale about thirty times smaller than the numerical box size. We therefore adopt an intermediate 
configuration with moderate value of the Reynolds number. 

It is not possible to evaluate the finite-size effects in the fully nonlinear case, especially because 
the theoretical study is based on additional assumptions. So we will restrict our analysis to the 
pure linear dynamics, or RDT. In so doing, we have to consider the following caveat: the pseudo- 
spectral method is assumed to be "exact" in the linear limit — to a given accuracy provided by the 
discretization in Fourier space — so that all RDT statistics derived from averaging u*Uj cannot be 
directly affected by the finite-size effect. On the other hand, statistics calculated from velocity 
components in physical space may be affected, even in the linear regime. 

In order to assess the influence of the confinement and the validity of DNS in the context of 
quasi-static MHD turbulence, we perform two simulations, for A'^ = 1 and 5, in which the nonlin- 
ear advective term is neglected. (Several comparisons of this type between linear and nonlinear 
simulations can be found in Favier et al. (2010).) These simulations can be compared to the linear 
analytical solutions from Rapid Distortion Theory (see Moffatt (1967) and Appendix B). In order 
to study specifically the effect of confinement, we compute correlation lengths defined by 



-{I) _ 
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{u.i{x)uj{x + r)) dr 



(13) 



where r^ = r5ki is the two-point velocity separation. In the current axisymmetric flow, the most 
relevant anisotropy indicators are the integral length scales with axial separation, relative to either 
axial or transverse velocity components (Cambon & Jacquin (1989)): 
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The expressions of linear solutions for these quantities can be found in Appendix B. As discussed 
above, these quantities evaluated by DNS are expected to coincide with these analytical formulas 
only in the theoretical limit of a projection base with an infinite number of degrees of freedom. 
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Figure 3: Results from EDQNM2 with different minimum wave number, from fcmjn = 0.01 to 
^min = 1- (a) Velocity correlation lengths. The horizontal dotted line corresponds to the numerical 
limit of 27r, present in DNS. (b) Angular energy spectra at t* = 12. 
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and L\-^^^ versus the dimensionless time t* = Iuq/Iq, 



Figure 2 shows the correlation lengths L33 
starting from the isotropic configuration at i = 0. The thick lines correspond to linearised DNS and 
thin lines correspond to linear analytical solutions. Both N ^ \ and iV = 5 cases are presented, 
for which one observes a growth of the correlation lengths, as expected in decaying turbulence. 
At moderate N , the length scales remain well below the numerical limit L = 27r, although one 
still notices a small departure between DNS and RDT results. The length scale predicted by 
analytical RDT is consistently larger than that of DNS. For N = 5, one clearly observes that the 
vertical correlation length L^j^ saturates before the maximum value 27r whereas the linear solution 
continues to grow. This difference is a clear example of confinement in anisotropic DNS and cannot 
be attributed to nonlinearities, which are absent in these simulations. 

Note that the initial integral length scale in the present DNS is very small {Iq w 0.25) compared 
to the computational box size. Removing completely all trace of numerical confinement would 
require decreasing Iq even more. Considering the current DNS resolution, the resulting Reynolds 
number would decrease too much for a turbulent flow to subsist. A solution is to increase the 
resolution, with increasingly demanding computational cost, to either a larger cubic box with 
resolution 2048^, or an adapted elongated box with resolution 512^ x 2048, as done by Vorobev et al. 
(2005) in MHD turbulence, or in rotating turbulence by Cambon et al. (1997) and in convective 
turbulence by Matsumoto (2009). The latter option indeed delays the confinement issue, which 
is most pregnant in the axial direction, but also implies to some degree the anticipation of the 
anisotropy in the later stage of the evolution. In the following, we retain a 512"^ resolution consistent 
with the isotropy of initial conditions, considering only the early time response t* < 6, hence 
keeping the flow in a significantly nonlinear regime while maintaining negligible confinement bias. 

In order to investigate whether the mechanisms observed at these low Reynolds numbers will 
persist at higher Rynolds numbers, we will use EDQNM closures. It therefore makes sense to 
address also the problem of confinement in the numerical resolution of EDQNM. Such confinement 
limitations should in principle also apply to EDQNM models since the minimum wave number is, 
as in DNS, fcmin = 1- However, the closure model is written in spectral space so that periodic 
boundary conditions are not explicit. The EDQNM spectral resolution can easily be increased in 
order to quantify the impact of numerical confinement through the value of the minimum resolved 
wave number kj^in- We thus perform three EDQNM2 simulations (the results are the same using 
EDQNMl) in the case TV = 5, with three different values fcmin = 0.01, 0.1, 1. Firstly, the time 
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evolution of the velocity correlation lengths L33 and L\-^ are plotted in figure 3(a). The predictions 
from the three simulations are almost undistinguishable, and, in contrast with the DNS results of 
figure 2(b), the growth of correlation lengths is not constrained by the value of the minimum wave 
number. Secondly, angular energy spectra, plotted on figure 3(b) at t* = 12, show that the spectral 
anisotropy is the same whatever kj^m (details on the anisotropic spectra will be presented in section 
6.2). Accordingly, we choose kmin = 1 in the following, to allow a complete comparison with DNS 
results, with the understanding that EDQNM is free from truncation effects. 
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Figure 4: Kinetic energy IC{t) and enstrophy < o;^ > /2 versus dimcnsionless time t* 
scale. — DNS, o EDQNMl, • EDQNM2, A isotropic EDQNM. 

6 Comparison between DNS and EDQNM 



in log-lin 



In this section, we propose a comparison between DNS and EDQNM results at moderate Reynolds 
number. In the context of quasi-static MHD turbulence, it is hardly possible with DNS to reach 
high Reynolds number simulations without encountering artificial efl^ects of the periodic boundary 
conditions, in view of the rapidly increasing numerical cost of pseudo-spectral methods with Re. In 
this section, we shall compare statistics obtained from the flow field predicted by DNS with data 
directly derived from EDQNM models, for a Reynolds number attainable by DNS. We first present 
a comparison of dynamical quantities in section 6.1, then an extended analysis of anisotropy in 
section section 6.2. 



6.1 Energetics 

Total kinetic energy and total enstrophy are presented in figures 4(a) and (b) respectively. After 
initialisation, the EDQNM model instantaneously builds triple correlations, or, in other terms, 
energy transfer spectra, close to the ones observed in DNS. Therefore, the initial dynamics for the 
EDQNMl-2 models and DNS are simflar. At larger times t* > 1, and for iV = 1, EDQNM2 {resp. 
EDQNMl) seems to overestimate {resp. underestimate) the value of kinetic energy. For N ~ 5, 
the kinetic energy and enstrophy decays predicted by EDQNM2 and DNS are in good agreement. 
In both cases, it appears that EDQNMl underestimates the kinetic energy and the enstrophy after 
the initial short time stage, whereas EDQNM2 predicts decay rates that are remarkably close to 
the DNS evolution, if one considers all the possible sources of statistical inaccuracies which may 
appear in DNS data. The good performance of EDQNM2 with respect to EDQNMl is clearly the 
sign that including the explicit effect of anisotropic Joule dissipation in the nonlinear dynamics is 
crucial for modelling quasi-static MHD turbulence. The scale-dependent Joule dissipation timescale 
Tuik) = I/A/q can thus be compared to the turbulent timescale r(fc) = £~^/^k~'^^^ , where £ is the 
kinetic energy dissipation. Equating these timescales yields a given wavenumber km = MqE"-^/^, 
say, which separates Joule dissipation dominated scales k < kM from dominant nonlinear dynamics 
k > kM (l/^Af is the equivalent of the Ozniidov scale introduced in stably stratified turbulence). 
For our runs at N — 1, fcjv/ — 5 initially, and kM — 50 at the end of the simulation, whereas for 
the run TV = 5, the figures are 60 at the beginning, and 800 in the end. This shows that, apart 
from the early stage of the iV = 1 case, in all our simulations, the energetic scales are dominated 
by ohmic dissipation (this is illustrated on figure 10). 

The axisymmetric EDQNM model is also valid for isotropic turbulence, but the numerical cost 
is considerably larger than that of the classical fully isotropic model. The results of isotropic DNS 
{i.e setting i?o = 0), presented on figure 4, are obtained from the same initial conditions, and 
show that the decay of kinetic energy is faster in the QS MHD case than in isotropic turbulence 
due to the additional ohmic dissipation. Concerning the evolution of enstrophy, one observes an 
initial increase for both EDQNM models in the isotropic case and in DNS, showing a short-time 
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Figure 5: Kinetic energy spectra at dinicnsionless times t* — 0, 3 and 6. 

re-adjustment which cannot occur when the interaction parameter N is large, since the magnetic 
effect catches up almost instantaneously. 

The kinetic energy spectra are plotted in figure 5 at three different times. The initial energy 
spectra are identical, since EDQNM spectra are initialized from DNS results. As already mentioned 
the initial integral length scale (resp. peak energy wave number) is smaller {resp. larger) than for 
classical hydrodynamic simulations. Figure 5 shows that the DNS and EDQNM2 spectral energy 
levels are in good agreement for all the scales of the flow. For iV = 1 (figure 5a), the slight 
overestimation of the energy by EDQNM2 is again observed, particularly at intermediate scales 
6 < A; < 12, while we retrieve the larger underprediction of the EDQNMl model. In all cases, the 
comparison between DNS and EDQNM in the dissipative range of the spectrum is not as good 
independently of the value of N and of the model version. Several explanations can be put forward, 
both on the account of the model or of the DNS approach: desaliasing in DNS, intermittency not 
present in the EDQNM model, truncation in both, etc. Overall, figure 5 still demonstrates that 
the EDQNM2 model is a good predictive model of the dynamics of QS MHD over a wide range of 
scales. 



6.2 Refined comparison of the anisotropy 

The level of anisotropy in the fiow can be quantified with increasing refinement degrees. A first 
measure is the ratio between horizontal and vertical kinetic energies, plotted in figure 6(a). The 
linear and inviscid regimes are characterised by the following scaling (Moffatt (1967)) 



2( 



u^) 



(16) 



where wy = Uz is the axial velocity component, and u± — \/y^ 
velocity components 



J-y, 



with Ur and m„ the transverse 



As already observed by Vorobev et al. (2005), Burattini et al. (2008a), 
Favier et al. (2010), this linear prediction is not observed in numerical simulations. The initial 

stage {t* < 1) is characterised by a decrease of the ratio r^ = (u^)/(un), in agreement with 

equation (16), but after a few turnover times, this ratio increases. It was shown that this is not 
due to a restoration of isotropy but to a nonlinear phenomenon linked to the particular quasi- 
two-dimensional structure of the fiow (for details, see Favier et al. (2010) and section 7.2). Figure 
6(a) shows that EDQNM2 reproduces this departure from the linear prediction, although with a 
time lag and a smaller amplitude. At small interaction parameter, EDQNM2 provides a better 
agreement with DNS than EDQNMl for iV = 1, less so for TV = 5. 

The ratio between transverse and axial kinetic energies presented on figure 6(a) sets the focus 
on the large scale dynamics. The small scale dynamics can be brought forward by computing a 
similar quantity based on vorticity components. We define the ratio between transverse and axial 
enstrophies as 
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Figure 6: (a) Ratio between horizontal and vertical kinetic energies versus dimcnsionlcss time 
t* — tuo/lo- (b) Ratio between horizontal and vertical enstrophies versus dimcnsionlcss time 
t* = tuo/la. — DNS, o EDQNMl, • EDQNM2, A Isotropic EDQNM. 

In a pure two-dimensional case, this ratio goes to zero, whereas in the isotropic case, it is about one. 
For A'' = 1 on figure 6(b), the ratio is always decreasing independently of the model considered, 
but is far from the two-dimensional value. For N = 5, there is a clear departure between DNS and 
EDQNM predictions. Initially in DNS, there is a strong decay of r^, then the trend is reversed 
synchronously with the decay reversal of r^ (figure 6a), at i* « 1. Eventually, r^^ decreases again. 
This three-stage evolution is not captured by the EDQNM model. The first increase stage after the 
initial decrease is reproduced, with a delay as for Tg, but the second change of slope is not. It seems 
that a phenomenon appears in DNS at i* w 2 — 3, whereby the ratio r„ decreases in DNS, which is 
not captured by the model. The multiplicity of possible nonlinear time scales in MHD turbulence 
might not be reproduced by the single time scale introduced in the closure (equation 33). 

Let us now compare the DNS and EDQNM results concerning the prediction of the directional 
anisotropy resulting from ohmic dissipation. The first effect of the magnetic field is to dissipate 
preferentially Fourier modes with wave vector k parallel to Bg. A direct consequence is the decrease 
of the transverse kinetic energy with respect to the axial one (since modes with k \\ Bq contribute 
only to transverse energy, see figure 1). This is observed in figure 6(a). 

The simplest way to quantify this directional anisotropy (directivity) is to consider typical 
angles defined in physical space, such as the ones introduced by Moreau and Shebalin (Alemany 
et al. (1979); Shebahn et al. (1983)). The 'Moreau angle' /3 defined by 



cos'^Pit)^ {JC{t)y 



cos^ ee{k,t)d^k, 



(18) 



directly derives from the one-point dynamical equation for the kinetic energy JC{t) — JJJ e(fc, i)d'^fe 

dIC/dt + 2Ml cos^ /3/C = -e (19) 

coming from integration of equation (9). This equation suggests as well to refine the definition of 
the separating wavenumber introduced in section 6.1 as /cm — M^ cos'^ /3e~^/^. 

The Shebalin angle, more widely used in the MHD community, characterizes the angular distri- 
bution of the vorticity spectrum k^e, as evidenced by its definition contrasted with equation (18): 



cos''0u{t)^{{io^)it))- 



k^ cos"^ 9e{k,t)d^k, 



(20) 



where the enstrophy is (w^) — JJJ k^e(k, t)d'^k. This definition is the continuous counterpart, in a 
slightly different form, of the classical discretized version Shebalin et al. (1983), used for the plots 
in figure 7(a): 

Efefcil«(fc,i)P 
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Figure 7: Anisotropic coefficients versus dimensionless time t 
(b) Anisotropic part of the Reynolds stress tensor 633 = b' 
EDQNM2. 
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+ bi3 



where k± — Jk"^ 
is the axial one. 



+ ky is the transverse component of the wave vector, and fcy = kz — k"^ cos^ 9 

On the other hand, we do not plot directly the Moreau angle here, but the 
equivalent quantity 6I3 defined by equation (22): 6I3 is proportional to the intensity of the first 
angular harmonic of e, through 5§3 = 1/6 — (1/2) cos^ (3. 

The Shebalin angles for the velocity field are first plotted in figure 7(a). In all cases — and 
similarly for the Moreau angles — , the increase from the isotropic initial value 6^ sa 54.7° in- 
dicates a concentration of energy in modes perpendicular to the imposed magnetic field. This 
two-dimensional limit corresponds to Ou ~ 90°. This is a well-known consequence of the ohmic 
dissipation, which results in physical space in a flow invariant in the axial direction. We note that 
EDQNM2 overpredicts the value of the Shebalin angle with respect to EDQNMl and DNS. This 
overestimation does not concern 633, as shown in figure 7(b). This suggests that the EDQNM pre- 
diction for the directional anisotropy is different for larger scales (energy distribution) and smaller 
scales (vorticity distribution), with a particular sensitivity of EDQNM2 at smaller scales. A small 
inaccuracy can therefore pull the Shebalin angle predicted by EDQNM2 in the wrong way, even if 
EDQNM2 gives a better overall prediction than EDQNMl. 

The polarization anisotropy is another kind of anisotropy that may appear in addition to the 
directional anisotropy. This anisotropy cannot be quantified with Shebalin angles, since it is not 
directly related to the dependence of the poloidal and toroidal velocity components upon 9, but is 
related to their difference. Its characterization requires a specific splitting of the deviatoric part 



6ij/^ of the Reynolds tensor Rij — {ui{x)uj{x)), where /C is the total kinetic 



bij = i?,j/(2/C) 

and 6ij the Kronecker tensor. Considering the axisymmetry of the flow about the axis of Bq, 

only one diagonal term is needed to describe the anisotropy, 633 say. Using equations (6) and (8), 



one obtains the two contributions for 633 
(1997)), with: 
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633 (Cambon & Jacquin (1989); Cambon et al. 
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^ f Z{k)sm^9d^k 
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where 9 is the polar angle between the wave vector k and the axis of symmetry (see figure 1), 
E(k) is the spherically-averaged kinetic energy spectrum, Z{k) is the polarization spectrum. As 
stated by its definition (22), ^33 is similar to the Shebalin angles in that it quantifies the directivity 
of the energy with respect to the vertical direction, bf^ quantifies the additional dimensionality 
anisotropy which is conveyed by the polarization spectrum Z. The limiting value b^^ = 1/6 is 
reached for two-dimensional flows, in both the 2D-3C and the 2D-2C cases, distinguished only by 
the value of b^: for 2D-3C flows, -1/2 for 2D-2C flows. 



Figure 7(b) presents the evolution of b 



33 



and 6|^3 



versus time. Concerning 633, the same con- 
clusions as the ones resulting from the Shebalin angles are drawn from the flgure. Note that it 
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Figure 8: Velocity correlation lengths for A^ = 1 (left) and A^ = 5 (right). Superscripts 3, h 
and i correspond to vertical direction (aligned with Bq), horizontal direction and any direction, 
respectively. — DNS, o EDQNMl and • EDQNM2. 




Figure 9: Evolution with time of C = ("D^is - 2{ul)L^^^. — DNS, o EDQNMl and • EDQNM2. 

is possible to rescale time with the ohmic dissipation characteristic time so that both 633 and 9^ 
collapse independently of the intensity of Bq (Favier et al. (2010) and analytical law in appendix B) . 

The polarization part 6^3, which is zero initially, decays in all cases, showing a global predomi- 
nance of toroidal over poloidal energy. Negative polarization is thus responsible for the increase of 
the componental enstrophy and velocity ratios r^ and r,. plotted in figure 6. EDQNMl-2 models 
underestimate the amplitude of polarization, but this is not necessarily a defect of the closure, given 
the spurious confinement effects yielding polarization in DNS, as shown and discussed in section 5. 
As already observed, EDQNM2 is in better agreement with DNS for A^ = 1 (when nonlinearities 
are important) whereas EDQNMl compares better for A" = 5 (when nonlinearities are dominated 
by ohmic dissipation). In view of the value of the separation scale kM presented in section 6.1, the 
dynamics is driven by nonlinear timescale only at the beginning of the simulation at A^ = 1. As 
mentioned in the Introduction, the additional physics injected into the EDQNM2 model through 
the straining timescale (see Appendix A) is corrected by the Joule dissipation time-scale, but the 
imbalance of the two acts variably depending on the regime. It seems here that the EDQNM2 
nonlinear improvements are too large for these Joule dissipation dominated scales. 

We then consider the time evolution of the velocity correlation lengths defined by equations (13), 

presented on figure 8. At the end of the simulations it* « 6), L33 « 2.2 and 0.98 for the respective 

(3) 

cases N = 5 and 1. The axial correlation length of axial velocity L33 is therefore always significantly 
smaller than the box size 27r. The axial correlation length of transverse velocity, however, for the 
case at A^ = 5, reaches about two thirds of the numerical box size. The correlation lengths obtained 
from EDQNM are close to the ones computed from DNS results, indicating a good prediction of 
the anisotropy of the large scale structures of the flow. As previously discussed in section 5, the 

(3) 

growth of L\^' computed by DNS seems to slow down in time, a fact that can be attributed to 



13 



the periodic boundary conditions. Such saturation is not apparent in EDQNM resuhs, so that 
the correlation lengths continue to grow. The confinement-related explanation is supported by the 
similitude between figure 2(b) and 8(b). 

Moreover, it is possible to isolate the contribution due to polarization by computing 



Jo 



kdk 



(24) 



This quantity is interesting for two reasons: (a) from equations (14) and (15), its departure from 
zero is only due to the polarization Z{k); (b) this quantity is accessible experimentally. Initially, 
C is exactly zero for EDQNM models since the polarization is set to zero at the beginning of the 
calculation. However, ({t* = 0) = —0.015 in DNS is small but not exactly zero. This may be 
a trace of the forcing scheme used to reach a quasi-steady state of hydrodynamic turbulence for 
t* < 0. C may also be dominated by contributions from small values of kz, where the DNS spectral 
discretization is too coarse to yield converged statistics. In all cases, figure 9 shows that C(^*) — C(0) 
decreases, in accordance with negative polarization. One observes that the equatorial polarization 
is underpredicted by EDQNM, which is consistent with the previous observations on the deviatoric 
tensor 633. However, the relative evolutions of the A^ = 1 and N — 5 EDQNM predictions for ( 
agree correctly with the dependence with N observed on the DNS curves. 

All the previous statistics involve a spectral integration over wave numbers, so that information 
about scale dependency is lost. On the contrary, the angular spectrum E{k, 6) retains both scalc- 
and angle-dependence: 
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where Afc and A6 specify the discretization steps in Fourier space used for computing the anisotropic 
spectra (see figure 1 in which the shaded region corresponds to the scales which contribute to 
E{k, 6)). Ring-averaged angular spectra E{k, 9) have already been used in the context of rotating 
turbulence by Cambon et al. (1997) and for stably stratified turbulence by Godeferd & Staquet 
(2003), and are similar to the ring decomposition by Burattini et al. (20086). We choose here 
Afc = 1 and IS.9 = tt/IO, figures that depend on the DNS resolution to ensure optimal statistical 
sampling. The angular spectra are plotted on figure 10, at time t* — 5. At the initial time t* — 0, 
all angular spectra collapse since the initial condition is isotropic. Figure 10 shows that, as time 
increases, most of the kinetic energy is concentrated in the spectrum with transverse wavevectors, 
since the Joule dissipation term in equation (9) reduces less energy at this orientation, indepen- 
dently on the wavenumber. The qualitative agreement of EDQNM model predictions with the 
DNS ones is impressive, considering the multi-scale, multi-directional character of these spectral 
statistics. There are, however, some differences. First, one observes that EDQNM2 overestimates 
slightly the equatorial kinetic energy, which is consistent with the overestimation of the Shebalin 
angle already observed in figure 7(a). However, the global angular dependency of the energy ob- 
served in DNS is well reproduced by EDQNM2, whereas EDQNMl overestimates the polar kinetic 
energy (see lowermost curves with o symbols on figure 10). In all models, as N increases, the 
angular anisotropy increases so that the flow tends to be invariant in the vertical direction. 

So far, we focused on the angular dependency of the kinetic energy. The departure from isotrop- 
ically distributed energy is due to Joule dissipation and is observable in the growth of the Shebalin 
angle 6'„ (see figure 7(a)), 6I3 (see figure 7(b)), and in angular spectra. However, it has been 
demonstrated that this effect is mostly linear, and that it can explain neither negative values of 
6^ (see figure 7(b)), nor the increase of the ratio between transverse and axial energies at large 
times (see figure 6). The poloidal/toroidal decomposition of spectral quantities (equation (6)), 
along with the angular spectral distribution, provides a way of understanding these unexplained 
features. Figure 11 presents the equatorial spectra (i.e. only transverse wave vectors are consid- 
ered) decomposed as poloidal (i.e. axial in this particular case, see figure 1) and toroidal {i.e. 
transverse in this configuration) contributions. The polarization anisotropy is clearly observable, 
as the difference between the two spectra. It is scale-dependent, with negative polarization at large 
scales ($^ > <I>^), responsible for the negative value of 6f3 and positive polarization at small scales 
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Figure 10: Angular energy spectra at i* = 5. ~ DNS, o EDQNMl and • EDQNM2. In each case, 
five curves are plotted, from top/equator {k _L Bq) to bottom/pole (fc || .Bo)- N = I (left) and 
A^ = 5 (right). 
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Figure 11: Equatorial energy spectra at i* — 6. The solid blue lines correspond to $(^)(fcj^), the 
dotted red lines correspond to $'^-'(/c_l). EDQNMl results are shifted down by three decades, 
EDQNM2 ones are shifted down by six decades. 

($1 < $2), The structure of the flow is therefore strongly scale-dependent with dominance of trans- 
verse kinetic energy at large scales and a dominance of axial kinetic energy at small scales. This 
departure from the poloidal/toroidal equipartition of energy is mainly observable for transverse 
wavevectors, where the energy accumulates because of ohmic dissipation. For axial wavevectors, Q 
goes to zero and this is no longer observable. Note that the cross-over wave number fc"]^ at which 
$^(/c^) = $^(fc^) (A:^ sa 20 on figure 11) depends mainly on the initial conditions and on the 
Reynolds number. Both EDQNMl-2 models reproduce this non-linear behaviour as well as the 
approximate location of the cross-over wave number. 

fc]] and fcj slopes are indicated on figure 11 for comparison with common scalings of two- 
dimensional turbulence with passive scalar (see Batchclor (1959), Bos et al. (2009), and the dis- 
cussion of the analogy with two-dimensional three components flows in section 7.2). 

6.3 Dynamical equilibrium and energy transfer spectra 

The anisotropic re-distribution of energy in quasi-static MHD turbulence, starting from isotropic 
initial turbulence, is the result of an essentially angular transfer, as we have shown above with DNS 
and the EDQNM model, and as was observed in towed-grid turbulence in mercury by Alemany et al. 
(1979); Caperan & Alemany (1985). These authors, using interaction parameters between N ~ 0.6 
and 1.17, also observe the appearance of a fcjT^ scaling for the axial kinetic energy spectrum i5||(fc||), 

— 5/3 

that progressively replaces the Kolmogorov scaling k., over an increasingly wider wavenumber 
range. The complete E\\{k±, fcii) distribution, plotted on figure 12, pictures the spectral equilibrium 
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Figure 12: Iso-contours of the two-dimensional spectral distribution of the axial kinetic energy 
£'ll {k±, fc||) in logarithm scale, (a) Figure extracted from Caperan & Alemany (1985), at Re ~ 1800 
and the interaction parameter N ~ 0.6. The conical spectral distribution is shown, (b) EDQNM2 
result at t* = 12, Re = 333 and A^ = 5, with iso-contour levels from -60 to -20 by steps of 5. 

of energy, due to both Joule dissipation — that drains energy towards the transverse 2D plane — and 
nonlinear inertial transfers. As argued by Caperan & Alemany (1985), the equilibrium between the 
two phenomena should lead to a conical distribution of spectral energy, which seems to be observed 
on figure 12(a). The same quantity computed with EDQNM is plotted on figure 12(b). The model 
permits this refined representation since it provides a smooth distribution of the spectra, hardly 
available in DNS. The comparison between the two panels of figure 12 suggests strong similarities 
in the dynamical equilibrium obtained in the experiment and in the EDQNM model. [From 
figure 10 which presents angular spectra, but contains the same information as shown differently 
on figure 12(b), we believe that an equivalent agreement would be obtained with DNS.] One must 
bear in mind, however, that the dimensional scalings of both plots of figure 12 are different, so 
that no quantitative agreement is claimed. 

In order to investigate further inertial transfers in the QS MHD turbulent flow, we compute 
energy transfer spectra. They are presented on figure 13 at the same time t* = 5 as the spectra 
of figure 10. In DNS, the spherically averaged transfer spectrum is directly computed from the 
nonlinear term s = u x u), with uj — 'V x u, as 
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where t ~ — fc^ [fc x (fe x s)]. We focus here on equatorial modes k _L Bq and we distinguish 
the axial equatorial transfer Ta{k±) and the transverse equatorial transfer Tt{k±). In EDQNM 
closures, these quantities are directly obtained as 



Tt (/c_l) - t'^"^ {k, 9 ^ tt/2) - t(^) (fc, 9 = 7r/2) 
Ta{kj_) = T^"^ (fc, 9 = 7r/2) + T^^^ (k, 9 ^ 7r/2) 



(27) 
(28) 



We observe an overall good agreement between DNS and EDQNM on figure 13. For A = 1, one 
observes a reduced transverse transfer compared to the axial one, both in DNS and in EDQNM 
closures. For N = 5, DNS and EDQNM2 clearly display a positive transfer at large scales, 
characteristic of an inverse cascade of kinetic energy. As described in Favier et al. (2010), the 
transverse component of the velocity behaves as in two-dimensional turbulence, with the axial 
velocity component acting as a passive scalar, thus characterised by a classical direct cascade. This 
inverse cascade of transverse velocity explains the reduction of dissipation and thus the dominance 
of transverse kinetic energy at large times (see figure 6). Note that EDQNMl is unable to reproduce 
the inverse cascade observed in DNS and EDQNM2. Finally, the oscillations observed in the DNS 
transfers for N = 5 could be explained by the fact that DNS yields one particular realization of the 
flow. The statistics of a flow can differ significantly from what is computed from an instantaneous 
flow field, in particular in the large scales. We therefore do not exclude that the double positive 
lobe of Tt{kj_) vanishes if we average over more flow realizations. 
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Figure 13: Equatorial kinetic energy transfer spectra at dimensionless times t* = h. — DNS, o 
EDQNMl and • EDQNM2. 

7 Additional results accessible only with the EDQNM clo- 
sure model 

We have presented in section 6 a comparison of the EDQNM2 closure model with DNS, which 
validates the results of the model for the given range of parameters attainable with DNS. However, 
due to the very way it is constructed and implemented, the added value of the EDQNM model 
is clearly to allow the investigation of an extended range of turbulent regimes. In the following 
two sections, we investigate high Reynolds number turbulence, currently out of the grasp of Direct 
Numerical Simulations (section 7.1), and a derived model for the limit case of two-dimensional 
three-components turbulence (section 7.2). 

7.1 High Reynolds number turbulence 

In this section, we address an important question of this article: are DNS predictions reliable to 
understand high Reynolds number quasi-static MHD turbulence given the moderate hydrodynamic 
Reynolds number? We use the EDQNM model at higher Reynolds number to answer this question. 
The number of wave numbers considered in EDQNM models has to be increased, along with the 
angular discretization and triadic interactions count. The following simulations are based on 100 
wave numbers, 48 polar angles and 48 angles for the direction of the plane of the triad around fe 
(denoted A in appendix A). The initial Reynolds number is increased from the previous value of 
Re ~ 333 up to Re w 2 x 10^. The initial condition for these high Reynolds simulations is similar to 
the one used for previous EDQNM simulations, except that the inertial range of the initial energy 
spectra is extended to higher wave numbers. 

Equatorial spectra are gathered on figure 14. First, the lowermost spectra on the figure recall 
the previous 512'^ DNS results (these spectra are shifted down by six decades). The corresponding 
cross-over wave number is fc^ ~ 171. The intermediate results correspond to EDQNM2, obtained 
with an initial Reynolds Re ~ 2200. The DNS resolution required to accurately simulate such a 
flow is about 1500^ Fourier modes. The results are qualitatively unchanged, but the slopes fcj^ 
and fcj] appear more clearly, even more so for the top curves on the figure, corresponding to an 
initial Reynolds number i?e ~ 2 x 10^. The corresponding DNS resolution using pseudo-spectral 
methods would be about 6000^ Fourier modes. We note also that the cross-over wave number fc^ 
defined by ^^{k'l) « ^^{k']_) increases with the value of the Reynolds number, to fc5_ — 35 for 
Re ~ 2200, and /c^ ~ 50 for i?e ~ 2 x 10^ 

7.2 A model for 2D-3C turbulence 

In two-dimensional three-components (2D-3C) flows, the velocity field contains three non zero 
components, which only vary in two directions (the transverse plane, say), and are independent of 



17 




le-16 



10 



100 



1000 



Figure 14: Poloidal/Toroidal decomposition of the equatorial energy spectra. The interaction 
parameter is A^ = 5. Ah the results are plotted at t* w 5. From bottom to top: DNS using 
512'^ Fourier modes (shifted down by 6 decades), EDQNM2 corresponding to a spectral resolution 
of 1500'^ Fourier modes (shifted down by 3 decades) and EDQNM2 corresponding to a spectral 
resolution of 6000^ Fourier modes. 
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Figure 15: (a) Ratio between horizontal and vertical kinetic energy for 2D-3C turbulence, (b) 
Anisotropic tensor 633 and its decomposition. — DNS and A EDQNM. 

the third direction (axial). The analogy between 2D-3C turbulence and the final state of quasi- 
static MHD turbulence is supported by theoretical (Montgomery & Turner (1982)) and numerical 
(Favier et al. (2010)) evidences. In previous sections, we found some indications, using EDQNM 
spectral closures, that this statement, supported by DNS at moderate Reynolds number, is valid 
for higher values of the Reynolds number, using EDQNM spectral closures. This last section is 
devoted to the comparison between DNS and EDQNM closures in a 2D-3C context. 

Theoretically, to consider 2D-3C turbulence is equivalent to considering purely 2D turbulence 
with a passive scalar (the latter being the vertical component of the velocity) . As shown by Cambon 
& Godeferd (1993) (see appendix A. 2), the EDQNMl model for anisotropic turbulence reduces 
exactly to a 2D-3C model for $^ and <i>^, in which $^ plays the same role as the scalar spectrum 
in 2D EDQNM (Lesieur & Herring (1985)). 

The previous 3D simulations tend to a 2D-3C state but this transition is triggered by dissipative 
effects so that the remaining energy is very small. To numerically investigate the 2D-3C state at 
high Reynolds numbers, we consider initially 2D-3C turbulence using both a 2D pseudo-spectral 
code and a 2D version of EDQNM closures presented above which include a passive scalar (consid- 
ered here as the axial velocity component). We use 1024^ Fourier modes for the DNS and 51 wave 
numbers for the spectral discretization of EDQNM. The initial condition is the same in both cases: 
^^{k,t* = 0) = <^'^{k,t* = 0) = 10-''fc2exp(-(fc/A:„)2), and fc,„ = 8. The molecular viscosity is 
fixed to i^ = 5 X 10~^ which corresponds to an initial Reynolds number of about 10'^. In the 3D 
axisymmetric case, the equatorial initial condition was also characterised by $"'^(A:_l) — ^^(fcj,), the 
main difference being that triple correlations were initially non zero. Here, the initial condition 
is a random Gaussian velocity field with an integral scale ^q ~ 0.32 and rms velocity mq ~ 0.18, 
hence with zero third-order moments. 

The ratio between transverse and axial kinetic energies is presented in figure 15(a). As expected, 
the initial value is about unity. As time increases, the inverse cascade of the horizontal velocity 
field develops so that the dissipation of horizontal components is reduced. This phenomenon is 
responsible for the growth of (m^ -I- Uy)/{ul). In the 3D axisymmetric case, one first observes a 
decrease of this quantity (see figure 6(a)). This is due to the transition from a 3D initial state 
to a quasi-two-dimensional state in which the inverse cascade occurs. Figure 15(a) also shows 
that the EDQNM evolution is faster (maybe from the fact that, in the EDQNM model, triple 
correlations — energy transfers — build up instantly). 

On figure 15(b), the anisotropic tensor 633 and its 633 , 633 decomposition are presented. Since 
the flow is 2D-3C, all axial derivatives are zero, so that all the energy is concentrated in the 
transverse plane. In that case, 6I3 reaches its maximum value 1/6 (see Gambon et al. (1997)). The 
polarization is initially very small, and becomes negative at larger times due to the dominance of 
toroidal (also transverse) energy with respect to the poloidal (also axial) energy. 

Finally, the poloidal/toroidal decomposition of the equatorial energy spectra is plotted in figure 
16. EDQNM and DNS are in very good agreement (again with the minor exception of the dissi- 
pative range). This result confirms the previous kj^ and kj^ scalings for the axial and transverse 
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Figure 16: Poloidal/Toroidal decomposition of the equatorial energy spectra for two and a half 
dimensional turbulence. All the results are plotted at t* « 5. — DNS, A EDQNM. 

velocity components, similar to those observed on figure 11 for the three-dimensional simulations. 
The asymptotic state of quasi-static MHD turbulence is therefore very similar to two-dimensional 
turbulence advecting a passive scalar. 

8 Conclusion and final remarks 

In this paper, we have investigated the dynamics and the detailed anisotropy of magneto-hydrodynamic 
turbulence in the quasi-static approximation at small magnetic Reynolds number, using Direct Nu- 
merical Simulations and a two-point statistical closure of EDQNM type. By essence, such closures 
consider statistical averages, which is a key advantage when considering turbulent flows, for two 
reasons: first, only one simulation is required to obtain averaged results, in constrast with the 
large number of realizations needed in DNS (typically more than a dozen) ; secondly, the obtained 
averages are smoother functions than in DNS, all the more if one considers high order moments {e.g 
third-order correlations). In terms of computational cost, isotropic EDQNM or the 2D-3C model 
presented in section 7.2 are thousands of times less costly than equivalent DNS. The axisymmetric 
anisotropic EDQNM2 model abandons one symmetry with respect to the isotropic context, thus 
the convolution integral is an order of magnitude more expensive. Therefore EDQNM2 computa- 
tions, although not as cpu and memory demanding as DNS by a factor of about 10 in the present 
parameter range, are also run on a parallel computer. The extension of two-point statistical clo- 
sures to bounded turbulent flows (Kraichnan (1972); Turner (2000); Laporta (1995)), however, is 
analytically and computationally challenging. 

In terms of statistical analysis, the closure allows for easy access to the general decomposition 
of tensors in the axisymmetric flow, such that refined statistics of turbulence can be used for 
characterizing anisotropy. The poloidal/toroidal decomposition of the velocity field and related 
second-order statistics permits the computation of a polarization tensor, which is a key indicator 
of whether the anisotropic mechanism is of linear nature — the Joule dissipation — or due to more 
complex nonlinear interactions. (The extraction of equivalent second-order statistics in physical 
space, although formally possible, would be hardly tractable, because of the differential operators 
involved). There remains the possibility to obtain such statistics by post-processing DNS data 
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fields, although with all the inaccuracies and sub-sampling issues due to limited resolutions. Clearly, 
DNS discretization is insufBcient in the very large scale range of the spectrum, and EDQNM is 
better off in this range and very adequate in the inertial range; less so in the smallest scales. 

We have nonetheless compared results of the EDQNM closure model with those of 512'^ DNS. 
The EDQNMl version of the model and the EDQNM2 one provide slightly different results, but 
the overall agreement with DNS is quite good. Comparisons involve kinetic energy and enstrophy, 
kinetic energy spectra and directional velocity correlation lengths. The latter allow us to address 
the question of numerical confinement due to the finiteness of the computational box in DNS, 
of importance in QS MHD turbulence in which the axial velocity correlation length increases 
tremendously. 

Several quantities were used to assess the level of anisotropy in the flow. Starting with initial 
conditions of isotropic turbulence, the ratios of transverse energy (resp. enstrophy) to axial energy 
(resp. enstrophy), the Shebalin angles and the off-diagonal components of the Reynolds stress 
tensor indicate all that the flow dynamics becomes closer to a two-dimensional three-components 
state. However, upon investigation of transverse and axial energy spectra, we are able to define a 
cross-over wavenumber below which the toroidal contribution dominates over the poloidal one, with 
a reversal of this order in the larger wave numbers or small scales. Not only are these predictions of 
DNS confirmed in a satisfactory quantitative manner by EDQNM, but the model allows to reach 
higher Reynolds numbers than permitted by DNS. The dynamics is not significantly altered at 
higher Reynolds numbers reached with the closure model. However, asymptotic scaling behaviour 
appears only very slowly. If a qualitative understanding of QS MHD is called for, both DNS 
and closure models are applicable. However, if scaling ranges and inertial range behaviour are of 
interest, two-point closures remain an indispensable tool. 

We conclude by noting that rotating turbulence bears strong similarities with QS MHD tur- 
bulence. In both cases, a transition from 3D to 2D structure is observed, and the 2D-2C trend is 
evidenced by the separation of L\^ and L33 integral scales, due to the growth of polarization in 
the horizontal transverse wave plane. This transition originates from the linear Joule dissipation 
term in QS MHD, but from nonlinear interactions dominated by cubic transfer terms such as T^'^' , 
when solid body rotation acts. Therefore, QS MHD turbulence may eventually become fully two- 
dimensional, whereas complete two-dimensionalization cannot be achieved in rotating turbulence 
in absence of additional phenomena. 

The authors thank the computing centre IDRIS of CNRS for the allocation of CPU time under 
project numbers 071433 and 022206. We also would like to thank the referees for they suggestions 
leading to improvements of the paper. 

A Detail on anisotropic EDQNM equations and their nu- 
merical calculation 

A.l EDQNM closure for the spectral energy transfers 

In section 3, the important term to specify is the quasi-normal one denoted 17^^,^,, (fc,p, t'), for 
modelling the fourth-order terms in (12) which is exactly given as a sum of quadratic terms from 
the set 

e = e(fc, t'),e' = e(p, t'),e" = e{q, t'),Z = Z{k, t') (29) 

in the case of a zero helicity flow. (The helicity, in contrast with the polarization anisotropy, 
remains zero if initially zero.) Instead of expressing ^g^'s''^ ^^ i^ simpler to derive its contribution 
to T^*^^ and T^^^ as was done for the EDQNM model for rotating turbulence, so that the numerical 
code for the EDQNM models used here is easily derived from the one for rotating turbulence (see 
e.g. Belief et al. (2006)). 

Detailed equations for T^^) and T*^) in the EDQNM2 model are 



2' s7^" ' ^ ^kpq + MS (cos2 9k + cos2 Op + cos2 e,) 
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Ai{sk, s'p, s"q)e"{e - e') + A2{sk, s'p, s"q)e^''"^" eZ{s"ci) + Aaisk, s'p, s"q)e^''^e" Z{sk) 

-A^isk, s'p, s"q)e'^''"^"e'Z{s"q) 
+A^isk, s'p, s"q) (e2'*"^"+2'^^Z(s"q)Z(sk) - e2'^"^"+2i^'^'z(s"q)Z(s'p))] d^p (30) 



and 

1 c r. / 



^•" = h^l 



s's" " ^ ^kpg + M^ (cos2 9k + C0S2 Op + cos2 0g) 



Asik, -s'p, -s"q)e"{e' - e) + yl4(fc, -s'p, -s" q)e'^''" ^" eZ{s" C{) + ^i(fc, -s'p, -s" q)e'^'^e" Z{]<i) 

-A^ik, -s'p, -s"q)e^'''^'e"Z{s'^) 
+A2{k, -s'p, -s"q) (e2'*"^"+2'-^Z(s"q)Z(k) - e^''" ^"+^'"'' ^' Z{s" qjZis'-p))] d^p (31) 



The geometric factors Ai to A^ are given in the appendix of Canibon et al. (1997), and in 
Sagaut & Cambon (2008); they depend only on the nioduh k, p, q, 'signed' by the polarization 
signs of helical modes, s — ±1, s' = ±l, s" — ±1. Ckpq depends only on the geometry of the triad 
as well, such that 

sin(p7g) _ sin(g, fc) _ sin(fc,p) _ 
k p q 

The internal triadic angles A, A' and A" denote the angle of rotation of the plane of the triad 
around k, p, q, respectively. Integration variables, which generate all the other terms at fixed fc, 
are p, q, as in isotropic EDQNM, and A, relevant in the axisymmetric case, and discretized as well. 
The only semi-empirical term in the formulae above is the viscous plus eddy damping term 
denoted 0^„„ since it is homogeneous to an inverse time scale, with 

%\ - iy{k^+P^ + q^)+^Hk,t)+^{p,t)+nq,t), (33) 

in which d{k,t) = A i J^^ p^E{p,t)dp] may be viewed as a straining decorrelation time scale 

of small turbulent structures by larger ones. Here, E{k) is the classical energy spectrum and 
A = 0.355 is the only adjusted constant of the model, computed from the Kolmogorov constant 
Ck with the relation Ck ^ 2.76^^/3 (Lesieur & Ossia (2000)). 

The EDQNMl version of the closure model does not incorporate linear Joule dissipation terms 
proportional to Mq in equations (30) and (31). It is therefore generic to any turbulent case, in 
which the distorsion only appears explicitly in a linear term added to the dissipation one. 

A. 2 Recovering the 2D-3C case 

This was done by Cambon & Godeferd (1993) as follows. In the 2D-3C limit, e and Z are concen- 
trated in the plane fcn (or k^ here) = 0, so that 

e(fc, t) == e(2^) (fc, t)5(fc|| ), Z{k, t) ^ Z(2^) (fc, i)(5(fc|| ), (34) 

and similarly for T^'^'^'. The Jacobian from {pi,P2) to {p,q) variables is now l/vT^^a?, only 
planar triads (fc|| = p|| = q^^ = 0) are called into play, and e^''^ — e^'^ — e^''^ = — 1. Accordingly, 
the 2D counterparts of Lin equation for $^ and $^ are derived as 



and 



{^^+2ukA<P\k,t)^T\k,t)^T^-'>^^''{k,t)-T^'''>-^''{k,t) 
fl + 2„kA $2(fc,t) ^ T2(fc,i) ^ T^-'>^^''{k,t) + T^''^'^''{k,t), 



(35) 
(36) 
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with 

and 

T\k,t)^ j j^ ^^;|(xy + z)$i(g,i)($^(p,t)-$2(fc,t))dpd<z. (38) 

It is shown that the 2D contribution from toroidal (horizontal in 2D) velocity is governed by the 
classical isotropic EDQNM equation restricted to 2D (Leith 1971, Pouquet et al. 1975), whereas 
the 2D contribution from poloidal (vertical in this limit) velocity is governed by the isotropic 
EDQNM equation in 2D for a passive scalar. 

More conventional relationship is found in term of the averaged spectrum using e^ '{k^t) = 
E{k,t)/{2Trk), as for the 3D isotropic case, in which e{k,t) = E{k,t)/{4:Trk'^). 

B RDT solutions for the correlation lengths 

The linear inviscid evolution of the spectral tensor is immediately found as 

e(fc,Ai,t) = ^^^^exp(-2MoVi), Zik,fi,t) = 0, (39) 

with /i — cos 9 and 9 the angle between k and the vertical. 

Two-dimensional energy components are invariant when defined as 

{ul){t)Li'^{t) = l/Co^o, {ul)it)L[\\t) =. {ul)it)Lil\t) = i/Co/o, (40) 

because they involve only contributions of e and Z at fx — 0. /Cq s-nd Iq are the initial kinetic energy 
and initial integral scale respectively. Kinetic energy and individual Reynolds stress components 
are given by 



/C(t) = /Co / exp (-2MoV^t) d^ 
Jo 



(41) 



and 



^0 /" '-■ ,2^„„„/ onj2,,,2\A., /„,2x/.N ''^0 



{4)it) = ^ y (1 - m') exp {~2M^t^,^) d^, {uj)it) = ^J (1 + ^2) ^^p (_2M2t^2) d^., (42) 

in agreement with d'^fe = 2Trk^dkd^ using polar-spherical coordinates for k and axisymmetry. 

The inviscid RDT time development of all relevant statistical quantities is derived analytically, 
in terms of the error function erf (exact relationship available from the authors upon request). 
The dominant terms in the evolution yield the following simple scalings: the kinetic energy decays 
as Mjj~^r((X))(2i)~^/^, as well as the Reynolds stress components; integral length scales with axial 
separation behave as MqIq^/I. Upon introduction of viscosity through the integrating factor e"^'''^ * 
in the integrands of equations (41) and (42), the viscous RDT solution is recovered, this time 
depending on the explicit shape of the spectrum E{k). For example, the RDT evolution of integral 
lengthscales may be compared to the evolution plotted on figure 2, and exhibit a linear evolution 
instead as the above inviscid ^/t behaviour. 

References 

Alemany, a., Moreau, R., Sulem, P. L & Frisch, U. 1979 Influence of an external magnetic 
field on homogeneous MHD turbulence. J. Mecan. 18, 277-313. 

Batchelor, G.K. 1959 Small-scale variation of convcctcd quantities like temperature in turbulent 
fluid. J. Fluid Mech. 5, 113. 

Bellet, F., Godeferd, F.S., Scott, J.F. & Cambon, C. 2006 Wave-turbulence in rapidly 
rotating flows. J. Fluid Mech. 562, 83-121. 

23 



Bos, W.J.T., Kadoch, B., Schneider, K. & Bertoglio, J. -P. 2009 Incrtial range scaling of 
the scalar flux spectrum in two-dimensional turbulence. Physics of Fluids 21 (If), ff5f05. 

BuRATTiNi, p., KiNET, M., Carati, D. & Knaepen, B. 2008a Anisotropy of velocity spectra 
in quasistatic magnetohydrodynamic turbulence. Phys. Fluids 20 (065ff0). 

BuRATTiNi, p., KiNET, M., Carati, D. & Knaepen, B. 20086 Spectral energetics of quasi-static 
MHD turbulence. Physica D 237, 2062-2066. 

Cambon, C. f990 Homogeneous MHD turbulence at weak magnetic Reynolds numbers: approach 
to angular-dependent spectra. In Advances in Turbulence Studies (ed. H. Branover & Y. Unger), 
Progress in Astronautics and Aeronautics, vol. 149, pp. 131-145. AIAA, Washington DC. 

Cambon, C. & Godeferd, F. S. f993 Inertial transfers in freely decaying rotating, stably- 
stratified, and MHD turbulence. In Progress in Turbulence Research (ed. H. Branover & 
Y. Unger), Progress in Astronautics and Aeronautics, vol. 162, pp. 150-168. AIAA, Washington 
DC. 

Cambon, C. & Jacquin, L. 1989 Spectral approach to non-isotropic turbulence subjected to 
rotation. J. Fluid Mech. 202, 295-317. 

Cambon, C, Mansour, N. N. & Godeferd, F. S. 1997 Energy transfer in rotating turbulence. 
J. Fluid Mech. 337, 303-332. 

Caperan, Ph. & Alemany, A. 1985 Turbulence homogenc MHD a faible nombre de reynolds 
magntique. etude de la transition vers la phase quasi bidimensionnellc et caracterisation dc son 
anisotropie. J. Meca. Theor. Appl. 4 (2), 175-200. 

Favier, B., Godeferd, F.S., Cambon, C. & Delache, A. 2010 On the two-dimensionalization 
of quasi-static MHD turbulence. Phys. Fluids 22, 075104. 

Godeferd, F. S. & Staquet, C. 2003 Statistical modelling and direct numerical simulations of 
decaying stably-stratified turbulence: Part 2: Large scales and small scales anisotropy. J. Fluid 
Mech. 486, 115-150. 

Herring, J.R. 1974 Approach of axisymmetric turbulence to isotropy. Phys. Fluids 17, 859-872. 

ISHIDA, T. & Kaneda, Y. 2007 Small-scale anisotropy in magnetohydrodynamic turbulence under 
a strong uniform magnetic field. Phys. Fluids 19 (075104), 10 pages. 

Jacquin, L., Leuchter, O., Cambon, C. & Mathieu, J. 1990 Homogeneous turbulence in 
the presence of rotation. J. Fluid Mech. 125, 505-534. 

VON Karman, T. & Lin, C. C. 1949 On the concept of similiarity in the theory of isotropic 
turbulence. Rev. Modern Phys. 21 (3), 516-519. 

Knaepen, B., Kassinos, S. & Carati, D. 2004 Magnetohydrodynamics turbulence at moderate 
Reynolds number. J. Fluid Mech. 513, 199-220. 

Knaepen, B. & Moreau, R. 2008 Magnetohydrodynamics turbulence at low magnetic Reynolds 
number. Annu. Rev. Fluid Mech. 40, 25-45. 

Kraichnan, R.H. 1972 Test-field model for inhomogeneous turbulence. J. Fluid Mech. 56, 287- 
304. 

Laporta, a. 1995 Spectral study and modelisation of an inhomogeneous turbulence. PhD thesis, 
cole Centrale de Lyon, in french. 

Lesieur, M. &: Herring, J. 1985 Diffusion of a passive scalar in two-dimensional turbulence. J. 
Fluid Mech. 161, 77-95. 

Lesieur, M. & Ossia, S. 2000 3d isotropic turbulence at very high Reynolds numbers: Edqnm 
study. J. of Turb. 1 (Art. No 7), 25 pages. 

24 



Matsumoto, T. 2009 Anomalous scaling of three-dimensional rayleigh-taylor turbulence. Phys. 
Rev. E 79, 055301. 

MOFFATT, H.K. 1967 On the suppression of turbulence by a uniform magnetic field. J. Fluid 
Mech. 28, 571-592. 

Montgomery, D. & Turner, L. 1982 Two-and-a-half-dimensional magnctohydrodynamic tur- 
bulence. Phys. Fluids 25(2), 345-349. 

Okamoto, N., Davidson, P. A. & Kaneda, Y. 2010 On the decay of Low-Magnetic- Reynolds- 
Number turbulence in an imposed magnetic field. J. Fluid Mech. 651, 295-318. 

Schumann, U. 1976 Numerical simulation of the transition from three- to two-dimensional tur- 
bulence under a uniform magnetic field. J. Fluid Mech. 74, 31-58. 

Shebalin, J. v., Matthaeus, W.H. & Montgomery, D. 1983 Anisotropy in MHD turbulence 
due to a mean magnetic field. J. Plasma Phys. 29, 525. 

Staplehurst, p. J., Davidson, P. A. & Dalziel, S. B. 2008 Structure formation in homoge- 
neous freely decaying rotating turbulence. J. Fluid Mech. 598, 81-105. 

Turner, L. 2000 Using helicity to characterize homogeneous and inhomogencous turbulent dy- 
namics. J. of Fluid Mech. 408, 205-238. 

Vorobev, a., Zikanov, O., Davidson, P. A. & Knaepen, B. 2005 Anisotropy of MHD tur- 
bulence at low magnetic Reynolds number. Phys. Fluids 17, 125105. 

Waleffe, F. 1992 The nature of triad interactions in homogeneous turbulence. Phys. Fluids A-4, 
350-363. 

Waleffe, F. 1993 Inertial transfers in the helical decomposition. Phys. Fluids A-5, 677-685. 

Zhou, Y. 2010 Renormalization group theory for fiuid and plasma turbulence. Phys. Rep. 488, 
1-49. 

Zhou, Y. & Matthaeus, W.H. 2005 Phenomenology treatment of magnctohydrodynamic tur- 
bulence with nonequipartition and anisotropy. Phys. Plasmas 12, 056503. 

Zhou, Y., Matthaeus, W.H. & Dmitruk, P. 2004 Magnctohydrodynamic turbulence and time 
scales in astrophysical and space plasmas. Rev. Mod. Phys. 76, 1015-1035. 



25 



